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ABSTRACT 

We present the results of a combined analysis of the kinematic and photometric proper- 
ties at large galactocentric radii of a sample of 14 low-luminosity early-type galaxies in 
the Fornax and Virgo clusters. From Gemini South GMOS long-slit spectroscopic data 
we measure radial profiles of the kinematic parameters v ro t, a, /13, and h± out to ~ 1 — 3 
effective radii. Multi-band imaging data from the HST/ACS are employed to evalu- 
ate surface brightness profiles and isophotal shape parameters of ellipticity, position 
angle and discyness/boxiness. The galaxies are found to host a cold and old stellar 
component which extend to the largest observed radii and that is the dominant source 
of their dynamical support. The prevalence of discy-shaped isophotes and the radial 
variation of their ellipticity are signatures of a gradual gas dissipation. An early star- 
forming collapse appears to be the main mechanism acting in the formation of these 
objects. Major mergers are unlikely to have occurred in these galaxies. We can not rule 
out a minor merging origin for these galaxies, but a comparison of our results with 
model predictions of different merger categories places some constraints on the possi- 
ble merger progenitors. These merger events are required to happen at high-redshift 
(i.e., z ^ 1), between progenitors of different mass ratio (at least 3:1) and containing 
a significant amount of gas (i.e., ^ 10 percent). A further scenario is that the low- 
luminosity galaxies were originally late-type galaxies, whose star formation has been 
truncated by removal of gas and subsequently the disc has been dynamically heated 
by high speed encounters in the cluster environment. 

Key words: galaxies: dwarf — galaxies: elliptical and lenticular, cD — galaxies: 
formation — galaxies: evolution — galaxies: photometry — galaxies: kinematics and 
dynamics 



1 INTRODUCTION 

In the past years, an increasing number of observa- 
tional studies have focused their effort in probing 
kinematic, photometric and stellar population proper- 
ties at large galactocentric rad ii of early-type gala x- 
ies 



is. 



Hau fc Forbes 20061 : ISpolaor et alT2009a; 



Proctor et al. 20091 

Foster et al. 2009). 



IWeiimans et al. 20091 : 



This interest is 



motivated by the 
fact that nuclear regions are indicative of only ~15 percent 
of the total stellar mass and of less than 10 percent of 
the angular momentum of a galaxy (i.e., assuming a de 
Vaucouleurs r 1 ' 4 light profile and a flat rotation profile). 

In numerical simulations, the efficiency of galaxy for- 
mation mechanisms are predicted to vary with radius leav- 
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ing measurable changes at different galactocentric radii. The 
spatial distribution of these observed stellar properties, and 
their relationship to galaxy structural parameters, provide 
unique clues to constrain competing galaxy formation pro- 
cesses. 

A specific ex ample is the mass-m etallicity gradient 
relation found by ISpolaor et al. f2009al ). The correlation 
shows a sharp change in slope at a dynamical mass of ~ 
3.5xl0 10 M (i.e., Mb ~ -19). The metallicity gradients of 
more massive galaxies (i.e., high- luminosity galaxies) flatten 
with increasing mass, perhaps due to the progressively larger 
number of mergers that these galaxies undergo to assem- 
ble their mass. In the remnants, the gradients are partially 
regenerated if a central burst of star formation is induced 
by the merger. In low- mass (i.e., low-luminosity) galaxies, 
an early star-forming collapse with a varying star forma- 
tion efficiency is required in order to explain the steepening 
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of gradients with increasing mass. The depth of the galaxy 
gravitational potential well coupled with energy feedback via 
supernova-driven galactic winds are responsible for shaping 
the gradients. 

The effect of gas dissipation, merger-induced star- 
bursts, and energy feedback are crucial ingredients 
in galaxy formation models to be able to reproduce 
the o bserved stellar proper ti es of early-type galaxies 
(e.g.. iKawata fc Gibson 20031; iPipino fc Matteucci 2006al ; 



ISpringel et al. 20051 ; iRobertson et al. 20061 ) . Specifically, 
the interplay and time ordering between dissipative 
processes and the redistribution of angular momentum 
could be the primary sources of the different stellar 
properties observe d in low-lumino sit y and high-luminosity 



galaxies fe.g.. iRix et al. 19991 ; Idj Matteo et al. 200 



I 



iHoffman et al. 2009^ 

The variety of kinematic and photometric features dis- 
played by early-type galaxies, such as isophotal shapes, 
extended stellar discs, amount of ordered and random 
motion are interpreted as signatures of these mecha- 
nisms. Fo r example, the r esults of the SAURON sur- 
vey (e.g., iBacon et al. 20011 ) have prompted the division 
of early-type galaxies into slow and fast rotators, de- 
pending on the amount of stellar a ngular momentum per 
unit mass inside one effective radius l|Emsellem et al. 20071 : 
iKrainovic et al. 20081 ). 

The aim of this paper is to compleme nt the stellar pop- 
ulations findings of ISpolaor et al. f2010bT ) (hereafter Paper 
II) on galaxy formation mechanisms acting in low-luminosity 
early-type galaxies. We do this by exploring the kinematic 
and photometric properties at large galactocentric radii of 
the same 14 low-luminosity galaxies. 

In general, properties of low-luminosity early- 
type galaxies have not been fu lly probed to large 
radii fe.g.. iDavies et al. 19831 : iBender et al. 19941 ; 



lHallidav et al. 20011 ). The low-surface brightness of these 
objects makes the measurement of reliable properties at 
large radii observationally challenging. Here, we study the 
galaxy kinematics by analysing the line-of-sight velocity 
distribution (LOSVD) as a function of galactocentric radius 
from high quality Gemini GMOS long-slit spectroscopic 
data. We derive spatially resolved radial profiles out to 
~ 1 — 3 effective radii for rotation velocity, velocity dis- 
persion and the Gauss-Hermite coefficients /13 and /14. The 
higher order moments of the distribution allow us to derive 
symmetric and asymmetric deviations from a pure Gaussian 
shape. They relate to the amo unt of ordered motion and the 
shape of stellar o rbi ts (e.g. van der Marel fc Franx 19931 ; 
IBender et al. 1994 ICappellari fc Emsellem 20041 ). The 
isophotal properties of a galaxy, such as surface brightness 
profile, ellipticity, position angle and discyness/boxiness are 
obtained from high-resolution HST/ACS imaging data in 
the g and z filters. 

The paper is organised as follows. In Section 2, we de- 
scribe the data sample. In Section 3, we describe observa- 
tions and reduction of spectroscopic and imaging data. Sec- 
tion 4 the technique used to recover stellar kinematics is 
described and the results are presented. In Section 5, we de- 
scribe the surface photometry analysis and the results are 
presented. In Section 6, the combined kinematic and pho- 
tometric properties of the galaxies are discussed. In Section 
7, we discuss our results in the context of predictions from 



competing galaxy formation scenarios. In Section 8, conclu- 
sions are given. 



2 DATA SAMPLE 

The data sample consists of 14 low-mass early-type galax- 
ies in the Fornax and Virgo clusters. The six galaxies 
in the Fornax cluster are chosen from the catalogue of 
iFerguson f 19891 ) . while the eight galaxies in the V irgo clus- 
ter are from the catalogue of lBinggeli et al. (19851 ). Galaxies 
of low mass are selected on the basis on their central stellar 
velocity dispersion ctq and the absolute B-band total magni- 
tude Mb , both of which are known to be independent proxies 
of galaxy mass. The sample uniformly covers the range 1.6 < 
log(ffo) < 2.15 and -16.5 > M B > -19.5. This translates to 
a dynamical mass range of about 10 9 < Mdyn/M© < 10 11 , 
using the log(Md yn ) = 21og(go ) + log(r e ) + 5.0 wher e r e is 
the effective radius in parsecs (ICappellari et al. 20061 ). The 
main properties of the galaxy sample are summarised in Ta- 
ble ED 



3 OBSERVATIONS AND DATA REDUCTION 

3.1 GMOS long-slit spectroscopy 

The observations were made using Gemini Multi-Object 
Spectrograph (GMOS) mounted on the Gemini South tele- 
scope in long-slit spectroscopy mode. The six galaxies in the 
Fornax cluster were observed during the semester 2006B 
(Program ID GS-2006B-Q-74), while the eight galaxies in 
the Virgo cluster were observed during the semester 2008A 
(Program ID GS-2008A-Q-3). 

Four 1800s exposure were taken for each galaxy ex- 
cept for FCC 335, for which two exposures of 1800s 
were taken. At the beginning of each exposure, the slit 
was centred on the nucleus and oriented along the major 
axis of the galaxy. The GMOS detector comprises three 
2048 x 4608 EEV CCDs with a pixel size of 13.5 x 13.5 )im 2 . 
We adopted the standard full frame readout, and to enhance 
the signal-to-noise, we selected to bin 2x spatially, yielding 
a spatial resolution of 0.1454 arcsec pixel . 

The spectrograph slit width was set to 1.0 arcsec for 
the 2006B run and 0.5 arcsec for the 2008A run. It was 
5.5' in length in both cases. For the 2006B run, we used 
the 600 line mm -1 grism blazed at 5075 A, providing a dis- 
persion of 0.457 A pixel -1 and a wavelength coverage of 
3670-6500 A. For the 2008A run, we used the 600 line mm -1 
grism blazed at 5250 A with 0.458 A pixel -1 , covering the 
wavelength range of 3860—6670 A. The instrumental set 
up provided a spectral resolution, as measured from the 
full width at half maximum (FWHM) of the arc lines, of 
~ 4.21 A (2006B run) and ~ 2.45 A (2008A run). Dur- 
ing the observations, the seeing varied from 0.7-arcsec to 
1.0-arcsec FWHM, which is more than adequate for the 
purposes of our study. Spectrophotometric and radial ve- 
locity standard stars covering a range of spectral classes 
were observed as calibra tors. These star s belong to both the 
Lick/IDS stellar library dWorthev 19941 ) and MILES library 
IjSanchez-Blazquez et al. 2 006). Bias frames, dome flat-fields 
and copper— argon (CuAr) arcs were also taken for calibra- 
tion. 
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Galaxy 


Alternative 


Hubble 


r e 


P.A. 


Distance 


M B 


M dyn 


Radial Range 




Name 


Type 


[arcsec] 


[degree] 


[Mpc] 


[mag] 


[10 9 M ] 


r e (kpc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


FCC 148 


NGC 1375 


SO 


14.9 


89 


18.45 


-17.30 


7.9 


2.0 (2.66) 


FCC 153 


IC 1963 


SO 


12.4 


81 


18.45 


-17.10 


3.1 


3.4 (4.02) 


FCC 170 


NGC 1381 


so 


13.7 


139 


18.45 


-18.88 


34.6 


3.4 (4.76) 


FCC 277 


NGC 1428 


E/S0 


9.8 


123 


18.45 


-17.52 


9.5 


2.0 (1.90) 


FCC 301 


ESO 358-G059 


E/S0 


10.4 


156 


18.45 


-17.02 


3.9 


1.0 (0.58) 


FCC 335 


ESO 359-G002 


E/S0 


15.6 


47 


18.45 


-16.93 


3.5 


1.0 (1.40) 


VCC 575 


NGC 4318 


E4 


7.2 


65 


22.08 


-17.61 


7.2 


2.1 (1.62) 


VCC 828 


NGC 4387 


E5 


9.1 


137 


17.94 


-18.26 


10.7 


2.0 (1.58) 


VCC 1025 


NGC 4434 


E0/S0 


12.9 


131 


22.39 


-18.72 


31.6 


2.0 (2.80) 


VCC 1146 


NGC 4458 


El 


26.1 


175 


16.37 


-18.13 


26.9 


0.8 (1.66) 


VCC 1178 


NGC 4464 


E3 


6.6 


11 


15.84 


-17.63 


7.2 


3.2 (1.62) 


VCC 1297 


NGC 4486B 


El 


2.3 


110 


16.29 


-16.77 


16.2 


3.0 (0.55) 


VCC 1475 


NGC 4515 


E2 


9.4 


15 


16.59 


-17.84 


7.7 


2.0 (1.51) 


VCC 1630 


NGC 4551 


E2 


12.5 


67 


16.14 


-18.06 


1.0 


2.2 (2.15) 



Table 1. Galaxy properties. Col.(l): galaxy name from the catalogue o nFerguson (19891) and iBinggeli et al. (19 85); (2): alternative 
galaxy name; (3): morphological type from the HyperLEDA database; (4): _B-band effective radius along the galaxy semi-major axis 
from the HyperLEDA database; (5): p osition angle of the major axis; (6): d istan ce to the galaxy using the distance modulus derived 
from surface bright ness fluctuations of iTonrv et al. (20011 1 , iMei et al. (2005T1 and iBlakeslee et al. (20091 ) ; we applied the correction of 
Uensen et al. (20031) to the values of Tonry; (7): total B-band absolute magnitudes estimated from total apparent magnitudes given in 
the NASA/IPAC Extra-galactic Database and corrected for Galactic extinctio n; (8): dynamical mass of the galaxy using the log(Md yn ) = 
2 log(cro) + log(r e ) + 5.0, where r e is the effective radius expressed in parsecs llCappellari et al. 2006); (9): radial range in effective radius 
and kpc (in parenthesis) covered by our observations. 



The data reduction was carried out using the Gem- 
ini/GMOS IRAF package (version 1.9.1). Each object frame 
is treated separately. Initial reduction of the CCD frames in- 
volves overscan correction, bias and dark subtraction, flat- 
field correction and cosmic ray removal. The Hatfield correc- 
tion is performed by means of both a quartz-halogen lamp 
(taken with the Gemini calibration unit, GCAL) and twi- 
light spectra, which are normalised and divided from all 
the spectra. This process corrects for pixel-to-pixel sensi- 
tivity variations and for two-dimensional large-scale illumi- 
nation patterns due to slit vignetting. Cosmic rays are iden- 
tified and eliminated by interpolating each line of an image 
with a high order function. The residuals (cosmic ray hits) 
that are narrower than the expect instrumental line width 
are replaced with the fitted value. The package GSWAVE- 
LENGTH is used to establish an accurate wavelength cal- 
ibration from copper— argon (CuAr) arc images taken just 
before and after each exposure. Each spectrum is converted 
to a linear wavelength scale using roughly 130 arc lines fitted 
by 3rd-5th order polynomials. An accuracy of better than 
0.1 A is consistently achieved for the wavelength calibration 
solution along the entire spectral coverage. The spectra are 
also corrected for geometrical distortions along the spectral 
directions (S-distortion). The spectral resolution is derived 
as the mean of the Gaussian FWHMs measured for a number 
of unblended arc-lamp lines which are distributed over the 
whole spectral range of a wavelength-calibrated spectrum. 

To obtain reliable measurements of the kinematic pa- 
rameters at large galactocentric radii, it is critical to perform 
an accurate sky subtraction. Sky contamination becomes in- 
creasingly significant for the outer parts of the galaxy. In 
these regions, the level of galaxy light can be a few per- 
cent of the sky signal. This effect represents one of the main 
sources of systematic errors in the analysis of the kinematic 



radial profiles. The generous field-of-view of GMOS and the 
intrinsic low luminosity of our sample allow us to accurately 
estimate the sky continuum from the edges of the slit, where 
the galaxy light is negligible. The sky level is estimated by 
interpolating along the outermost 10—20 arcsec at the two 
edges of the slit and subtracted from the whole image. A 
sky subtraction accurate to within one percent is achieved, 
guaranteeing negligible related systematics on the measured 
kinematic parameters. 

For each galaxy, we obtained four fully reduced galaxy 
frames of identical exposure time. These are then co-added 
to form a single frame. The final galaxy frame is a two- 
dimensional spectrum from which we extract ID spectra 
along the slit. The spatial width (i.e. the number of CCD 
rows binned) for each extracted ID spectrum increases with 
radius to achieve a minimum signal-to-noise (S/N) per A of 
~ 35 in the spectral region of Mgb. By adopting this crite- 
rion we have been able to obtain reliable measurements of 
kinematic parameters to galactocentric radii of ~ 3r e . 



3.2 HST/ACS imaging 

The imaging observations were performed with the Ad- 
vanced Camera for Surveys (ACS) mounted on the Hub- 
ble Space Telescope (HST). The eight Virgo cluster galaxies 
were observed by the ACS Virgo Cluster Survey (ACSVCS; 
ICote et al. 20041 ). while the six galaxies in the Fornax cluster 
were observed by t he ACS Fornax Cluster Survey (ACSFCS; 
I Jordan et al. 20071 ). Each galaxy was ob served using the 
ACS Wide Field C hannel (WFC) mode (|Ford et al. 19981 ; 
ISirianni et al. 20051 ). The WFC consists of two 2048x4096 
CCD detectors (WFC1 and WFC2) with a pixel size of 15^m 
and a spatial scale of 0.049" pixel -1 . The nucleus of each 
galaxy was placed on the WFC1 detector at the reference 
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pixel (2096, 200). To avoid the gap between the WFC1 and 
WFC2 detectors an additional offset 5" to 10" was applied 
(depending by the brightness of the target galaxy) . 

Each galaxy was observed with the same observational 
methodology. Two images in the F475W filter (which resem- 
bles the SDSS g filter) for a total of 760s, with a respective 
exposure time of 380s to aid the removal of cosmic rays. 
Two images in the F850LP filter (similar to the SDSS z fil- 
ter) for a total of 1220s, split in a pair of 760s exposures. 
A further image of 90s exposure in F850LP filter was taken 
to recover the central regions of the galaxies whose nuclei 
might saturate in the longer 560s exposure. 

The data were downloaded from the Multimission 
Archive at the Space Telescope Science Institute (MAST). 
The image s were processed by th e ACS calibration pipeline 
(GALACS; |Pavlovskv et al. 20061 ') using the "On the Fly Re- 
processing" option. An automated dat a reduction pipelin e 
was developed for the ACSVCS data l| Jordan et al. 2004h . 
The data reductions procedures for the ACSFCS data have 
been documented in I Jordan et al. (20071 ) and they are simi- 
lar to those adopted for the ACSVCS data. Briefly, CALACS 
performed bias and dark subtraction, flat-fielding, and over- 
scan removal. The five images of each galaxy were aligned by 
matching coordinates for a list of objects in common, and 
images taken with the same filter were then combined. A 
cosmic rays rejection procedure was applied before perform- 
ing geometric distortion correction to combined images. The 
field-of-view of the fully reduced images is approximately 
202" x 202" in the shape of a rhomboid. 



4 GALAXY STELLAR KINEMATICS 

We use the penalised pix el fitting (pPXF) code of 
ICappellari fc Emsellem (20041 ) to extract kinematic param- 
eters of the galaxies. This is based on the method for iden- 
tifi cation of non-Gaussian line pr ofiles, originally developed 
by Ivan der Marel fc Franx f 19931 ). 

In the pPXF software the line-of-sight velocity distribu- 
tion (LOSVD) is parametrized as a sum of orthogonal func- 
tions in a Gauss-Hermite series, namely a Gaussian plus a 
series of Gauss-Hermite polynomials. A \ 2 minimising pro- 
cess determines the best-fitting moments of the LOSVD. 
The mean line-of-sight velocity, v\ os , the velocity dispersion, 
a, and the coefficients of the Gauss-Hermite polynomials of 
third, /13, and forth, /14, order. 

The coefficients /13 and /14 quantify asymmetric (skew- 
ness) and symmetric (kurtosis) deviations of the LOSVD 
from a simple Gaussian distribution, respectively. Positive 
/13 values indicate a distribution skewed towards velocities 
lower than the systemic velocity of the galaxy, v 3ys - Con- 
versely, negative values of h% are characteristic of a distri- 
bution skewed toward velocities greater than v sy3 . Positive 
hi values correspond to a distribution more peaked than a 
Gaussian, while negative values of hi denote a flatter distri- 
bution. 

The software works in pixel space, finding the linear 
combination of template stars spectra which, after being 
convolved with an appropriate LOSVD, best reproduces the 
galaxy spectrum. The advantage of this technique is that 
errors due to template-mismatch are reduced by using a 
combination of template stars. In our analysis, ten template 



stars of spectral type in the range from G3 to K9 are used. 
The entire wavelength range covered by the galaxy spectrum 
is modelled, but the code allows us to mask out regions pos- 
sible affected by skyline residuals or emission lines. A pixel- 
by-pixel x 2 minimisation of the residuals to the best fit is 
carried out. An adjustable penalty term is applied to the 
X 2 to bias the higher order moments coefficients h$ and /14 
of the resulting LOSVD towards a Gaussian shape. We find 
that a small bias of 0.4 provides an optimum value for the 
spectra of our data sample. 

We estimate random errors on the derived parameters 
using a series of Monte Carlo simulations. The entire kine- 
matic measurement process is repeated on a large number 
of different realisations of our data sample. The model spec- 
tra are obtained using the corresponding best linear com- 
bination of template stars, broadened by the best-fitting 
LOSVDs and degraded to a range of S/N values. 



4.1 Radial profiles 

The radial profiles of the kinematic parameters v Iot , u, h-j, 
and /14 for the 14 low-mass galaxies are presented in Fig. [l] 
Fig. [2] and Fig. [3] The parameters are expressed as a func- 
tion of the galactocentric radius along the galaxy's major- 
axis of each spatially resolved aperture, scaled with the ef- 
fective radius of the galaxy's semi-major axis. The rotation 
velocity v Iot is calculated as the difference between the mea- 
sured v\ os of an aperture and the v sys of the galaxy, that is 
the median value of v\ OB from all apertures extracted along 
the slit for each galaxy. A description of the individual galax- 
ies is given in Appendix |A1 

To be able to compare our results with previous studies 
we derive the following parameters describing the kinematic 
properties of the galaxies. The maximum measured rota- 
tion velocity w ma x within a radius r ^C r e /2. The central 00 
and mean a velocity dispersion, obtained by averaging the 
a values estimated in all spatially resolved apertures within 
a radius r ^ r e /8 and r ^ r e /2, respectively. The mean of 
absolute values of /13, | /13 j, and the mean of values of /14, 
/14, estimated from all spatially resolved apertures within a 
radius r ^ r e /2. The parameters for each galaxy are listed 
in Table [2] 



5 GALAXY SURFACE PHOTOMETRY 

The isophotal analysis of images in both filters of each 
galaxy is conducted using the I RAF task ELLIPSE . The 
package is based on the method of lJedrzeiewski (19871 ). and 
it involves the creation of a smooth elliptical galaxy mod el 
to fit to the real galaxy image fe.g.. lFerrarese et al. 20061 ). 

The technique makes the assumption that each isophote 
can be modelled by an ellipse whose centre, ellipticity and 
position angle are allowed to vary to reach the best-fitting 
values. The model is built up from a first guess elliptical 
isophote whose centre is pre-determined. The semi-major 
axis of each isophote is then incremented logarithmically 
from the previously fitted ellipse. The process is iterated 
until specific convergency criteria are achieved. A standard 
k-sigma clipping routine is used to evaluate and to flag diver- 
gent points, thus increasing the stability of the convergency 
process. A mask is created to exclude from the fit bright 
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Figure 1. The major-axis stellar kinematics of the Fornax cluster low-luminosity galaxies. From bottom to top: radial profiles of 
rotational velocity v ro t, velocity dispersion a and the Gauss-Hcrmitc coefficients /13 and /14. The kinematic quantities are expressed as a 
function of the galactocentric radius of each spatially resolved aperture, normalised with the effective radius of the galaxy's semi-major 
axis. 
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Figure 2. The major-axis stellar kinematics of the Virgo cluster low-luminosity galaxies. From bottom to top: radial profiles of rotational 
velocity v ro t, velocity dispersion a and the Gauss-Hermite coefficients hz and /14. The kinematic quantities are expressed as a function 
of the galactocentric radius of each spatially resolved aperture, normalised with the effective radius of the galaxy's semi-major axis. 
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Figure 3. The major-axis stellar kinematics of the Virgo cluster low-luminosity galaxies. Continued from Fig. [2] 
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Table 2. Kinematic and photometric parameters for the 14 low-luminosity early-type galaxies. See text for details. 



globular clusters, foreground stars, background galaxies and 
the gap between the individual CCD detectors. 

The background-subtraction is based on the average of 
median values from several 10x10 pixel boxes sampled in 
regions not affected by galaxy light. The conversion from 
the HST/ACS instrumental F475W and F850W filters sys- 
tem to the AB photometri c system is performe d using the 
synthetic coefficients from ISirianni et al. ('20051 ') . In partic- 
ular, the F475W and F850W filters are transformed into 
qab and zab , respectively. Final magnitudes were corrected 
for Galactic extinctio n using the reddenin g values from the 
DIRBE dust maps of lSchlegel et al. (19981 ). 

The errors in magnitude are computed from the rms 
scatter of intensity data along the fitted ellipse and the er- 
rors in the background subtraction. Errors in the Galactic 
extinction correction values and in the zero-points are neg- 
ligible with respect to the errors in the intensity profiles. 
Ellipse geometry parameter errors are obtained from the in- 
ternal errors in the harmonic fit. Harmonic amplitude errors 



are obtained from the fit errors after removal of all harmon- 
ics up to and including, the one being considered. 



5.1 Radial profiles 

We present the results of the isophotal analysis for the 14 
sample galaxies in Fig. [4] Fig. [5] Fig. [6] and Fig. [7] Pho- 
tometric parameters are plotted as function of the geomet- 
ric mean radius r, defined as r — a[l — e(a)] 1 ' 2 where a 
is the length of the isophotal semi-major axis. This radius 
is then scaled with the galaxy's effective radius along the 
semi-major axis and expressed as r/r e . 

For each galaxy, the /i 9 , /i z surface brightness profiles 
and the g - z colour index profile are shown. We also plot 
radial profiles of isophotal geometry parameters of elliptic- 
ity e, position angle P. A. and the isophote's deviations from 
perfect ellipticity A3, A±, and B3, B4. These parameters are 
respectively the third and fourth cosine and sine terms of a 
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Figure 4. Isophotal properties of the Fornax cluster low-luminosity galaxies. Open black squares indicate g-band parameters. Filled 
black squares denote z-band parameters. In each plot are shown: surface brightness profiles fj, and colour index profiles; ellipticity e, 
position angle P. A. and best fit parameters A3, A4, and S3, B4. Positive values of B4 indicate discy isophotes, whereas negative values 
represent boxy isophotes. The photometric quantities are expressed as a function of the galactocentric radius of data points, normalised 
with the effective radius of the galaxy's semi-major axis. 



Fourier expansion used by the software. The coefficients A3 
and -B3 indicate isophotes with three-fold deviations from 
ellipses, i.e. egg-shaped or heart-shaped. Four-fold devia- 
tions are represented by the A4, and B4, coefficients, such 
that a rhomboidal or diamond shaped isophote translates 
into a nonzero A4, value. A positive value for B4 indicates 
discy isophotes, whereas a negative value represents boxy 
isophotes. 



values measured inside half the effective radius in both fil- 



ters. Mean parameter values of colour index g — z, ellipticity 
e, position angle P. A., coefficient B4 for all 14 galaxies are 
listed in Table 



Photometric parameters are derived by averaging all 
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Figure 5. Isophotal properties of the Fornax cluster low-luminosity galaxies. Continued from Fig. [4] 



6 GALAXY PROPERTIES 

6.1 Vrot/u radial profiles 

The radial profiles of the ratio between the local rotation ve- 
locity and local velocity dispersion for our 14 low-luminosity 
galaxies are shown in Fig. [SJ and Fig. [5] This ratio allows 
us to assess the level of ordered versus random motion as 
a function of galactocentric radius. The profiles have been 
folded onto the positive side of the major axis, with a sign 
change for the rotation velocity. 

For each galaxy, we find that the profiles of the two 
major axis's sides are similar within the errors. This suggests 
an overall symmetry between the two sides, and thus a large- 
scale kinematic coherence. 

In general, the profiles raise with the distance to centre 
of galaxy. Inside one effective radius the value of v ro t/a ratio 
covers the range - 1.4. At two effective radii the galaxies 
have v/a ratio values in the range 0.6 - 2.4. This is attributed 
to the rotation velocity increasing with radius and contem- 
porarily the velocity dispersion decreasing. In other words, 
the galaxy kinematics is primarily driven by the ordered 
motion of stars. 



6.2 Rotational support versus Anisotropy 

We present the v max /a — e diagram IjBinnev 197a ) for our 
14 galaxies in Fig. 1101 The plot allows us to quantify the 
amount of rotational support of these objects with respect 
to their intrinsic shape. We adopt as variables the ratio of 
ti mal to the mean stellar velocity dispersion a and the mean 
isophotal ellipticity e. 

In the diagram we also sho w the relation ship for oblate 
isotropic rotator models from iBinnev (19781 ). The model 
predictions are derived from the tensor virial theorem, as- 
suming that elliptical galaxies are oblate spheroids of con- 
stant ellipticity and with isotropic velocity distributions. 



The values measured for our galaxies are generally consis- 
tent with the theoretical relationship, thus being oblate ro- 
tationally flattened objects. In the plot we also show the 
data for low-luminosity elliptical galaxies in t he Virgo clus- 
ter and Leo group from lHallidav et al. (20011 ). Four of their 
14 galaxies are in common with our data sample. Their mea- 
surements are in good agreement with ours, however, we 
decide to use our values. 

In Fig. 1111 we plot the anisotropy parameter (v/a)* 
against the total B-band absolute magnitude for each 
galaxy. The anisotropy parameter expresses the ratio 
between the observed value of Vmax/a and the (v/a)oi 
value expected for an oblate rotator (jBinnev 197a ; 
iDavies et al. 19831 ). We define (v/a)* in terms of observable 
quantities as (vmax/ a)/ \/t/(l — I) (e.g. IfCormendv 19821 ; 
iBender fc Nieto 1990t ). We also pl otted the data fo r 
dwarf elliptical ga laxies from IPedraz et al. ( 2002). 



iGeha et al. (2003[) and van Zee et al. (20041) . low- luminosity 
ellipticals in lHallidav et al. (200 ll) and h igh-luminosity el- 
liptical galaxies from Be nder et al. (19941 ). 

The anisotropy parameter is independent of ellip- 
ticity and projection effects, thus it can be used to 
separate galaxies with different dynamic al support (e.g., 
IDavies et al. 19831 : IScorza fc Bender 19951 ). Objects with 
(v/a)* < 0.7 are anisotropically supported systems, whereas 
values of (v/a)* ^ 0.7 define objects primarily sup- 
ported by rotation. Despite the scatter, the results suggest 
that high-luminosity elliptical galaxies are predominantly 
anisotropic objects and the amount of rotational support 
becomes increasingly significant towards intermediate and 
low-luminosity galaxies. The rotational component of the 
latter is significantly strong such that their observed flat- 
tened morphology is attributed to fast rotation. Dwarf el- 
liptical galaxies tend to deviate from this trend, although 
exhibiting a range of rotational properties. 
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Figure 6. Isophotal properties of the Virgo cluster low-luminosity galaxies. Open black squares indicate g-band parameters. Filled black 
squares denote 2-band parameters. In each plot are shown: surface brightness profiles fi and colour index profiles; ellipticity e, position 
angle P. A. and best fit parameters A3, A4, and B3, B4. Positive values of B4 indicate discy isophotes, whereas negative values represent 
boxy isophotes. The photometric quantities are expressed as a function of the galactocentric radius of data points, normalised with the 
effective radius of the galaxy's semi-major axis. 



6.3 Rotational support versus Discyness/Boxiness 

In Fig. 1121 the relation between amount of rotational sup- 
port and isophotal shape is examined. We plot anisotropy 
parameter against the B4 coefficient of our sample galaxies , 
and the values of bright ellipticals by iBender et al. (1994 ) 
for comparison. 

There exists a relation such that galaxies with discy 
isophotes tend to be more rotationally supported than boxy 



galax ies. The relation is independent of galaxy's luminosity 
(e.g., IBender et al. 1994 ). We find that almost 80 percent 
of the low-luminosity galaxies have discy-shaped isophotes 
or in other terms an excess of light along the galaxy's ma- 
jor and/or minor axis. The remaining objects have boxy 
isophotes, for which the excess of light lies along a line at 
45° with respect to the galaxy's axes. 
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Figure 7. Isophotal properties of the Virgo cluster low-luminosity galaxies. Continued from Fig. EJ 



6.4 LOSVD deviations versus Rotational support 

The plot of the skewness coefficient hs against the rota- 
tional support parameter v rot /<J is shown in Fig. 1131 Each 
point corresponds to a local position along the galactocen- 
tric radius, and data for all 14 low- luminosity galaxies are 
presented in the same diagram. 

The two parameters define a negat ive correlation, ofte n 
referred as the h-j, — v rot /a plane (e.g. iBender et al. 1994T) . 
In particular, in the central regions of galaxies, which cor- 
responds to low \v ro t/cr\ values, the coefficient | /zs | increases 
rapidly with \v ro t/o'\- At high values of \v ro t/cr\, characteris- 
tic of intermediate and large radii, \hz\ reaches its maximum 



values, up to \hz\ ~ 0.15. The relation was firstly observed by 
IBender et al. fl994T) . which noted that its slope is too steep 
to be explained by two- in tegral oblate rotators models (e.g., 
iDehnen fc Gerhard 1994 ). Instead, they suggested that el- 
lipticals harbouring large-scale stellar discs (~ 15 percent 
in mass, scalength ~ r e ) could explain the I13 — v/a plane. 
Our results for low-luminosity galaxies confirm the validity 
at large radii r > lr e of the hz — v ro t/<J distribution observed 
by IBender et al. (1991 ). 

We now turn to the kurtosis parameter hi, estimat- 
ing the symmetric deviations of the LOSVDs from a Gaus- 
sian shape. In Fig. [14] we present the hi — v rot /<J dis- 
tribution for all the low-luminosity galaxies. In contrast 
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Figure 8. Radial profiles of the rotational support parameter 
(ratio of local rotation velocity to local velocity dispersion) for 
the six Fornax cluster galaxies. The quantity is expressed as a 
function of the ratio between the local radius and the galaxy's 
effective radius r/r e . Continuous and dashed lines indicate the 
two sides with respect to the galaxy centre. 



Figure 10. Ratio of maximum rotation velocity v m ax to mean 
velocity dispersion <r plotted against mean isophotal ellipticity e. 
Filled red squares are the data for the 14 low-luminosity galax- 
ies presented in this paper. Open black squares are values by 
lHallidav et al. (2001! ) for low-luminosity elliptical galaxies. The 
solid black line represents the expected relationshi p for an oblate , 
isotropic model galaxy flattened by rotation from lBinnev ( 1978). 
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Figure 9. Radial profiles of the rotational support parameter 
(ratio of local rotation velocity to local velocity dispersion) for 
the eight Virgo cluster galaxies. The quantity is expressed as a 
function of the ratio between the local radius and the galaxy's 
effective radius r/r e . Continuous and dashed lines indicate the 
two sides with respect to the galaxy centre. 
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Figure 11. Logarithm of the anisotropy parameter (v/cr)* (ra- 
tio of the observed value of v m ax/cF to the value predicted 
by the isotropic model) plotted against the total B-band ab- 
solute magnitude. Filled red squares are the data for the 14 
low-luminosity galaxie s presented in this paper. Open black 
squares are values by lHallidav et al. (20011 ) for low-luminosity 
elliptical galaxies. Open black circles, open black triangles and 
filled black trian gles are values for d warf elliptical ga laxies 
respectively from I van Zee et al. (2004r ). iGeha et al. (20031 ) and 
IPedraz et al. (20021). Filled black circle s are values for bright el- 
liptical galaxies bv lBender et al. (19941 ). 
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Figure 12. Anisotropy parameter against the B4 coefficient. 
Filled red squares are the data for the 14 low-luminosity galax- 
ies presented in th is paper. Open black circles are values by 
iBender et al. (19941 ) for bright elliptical galaxies. Positive values 
of B4 indicate discy isophotes, whereas negative values represent 
boxy isophotes. 
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Figure 13. Local relation between the Gauss-Hermite coefficient 
/i3 and the rotational parameter v/a. Data points are from the 
14 low-luminosity galaxies presented in this paper. 



to /13, /14 is not observed to correlate with the rotational 
parameter. This lac k of correlation was already found by 
IBender et al. f 19941 ). However, the distribution present a 
wide spread of positive ft.4 values at any v rot /o, indicating 
a significant rotationally-supported disc component. 



0.2 



0.1 



-0.1 



* * 



-2-10 1 2 

V rotA 

Figure 14. Local relation between the Gauss-Hermite coefficient 
/14 and the rotational parameter v/a. Data points are from the 
14 low-luminosity galaxies presented in this paper. 



7 DISCUSSION 

We find that almost 85 per cent of our sample galaxies are 
characterised by a cold stellar component extending out to 
the outermost regions reached by our study, i.e. ~ 3r e . The 
stellar velocity dispersion is observed to quickly decline with 
radius, such that ordered motion prevails over random stel- 
lar orbits. The galaxies are well described by oblate rotation- 
ally supported models. The results of our analysis confirm s 
previous circumsta ntial fi ndings by IBender fc Nieto (199(1 ). 
iDavies et al. (19831 ), and lHallidav et al. (200 j ). However, 
unlike these earlier studies, we are able to spatially re- 
solved enough apertures at large galactocentric radii (i.e., 
r ^ 0.5r e ), where the data have a sufficiently high S/N to 
make robust measurements. 

During an early star-forming col lapse (e.g., 
IChiosi fc Carraro 20021 ; iKawata fc Gibson 20 03) , a cold and 
extended stellar component can naturally originate from the 
high angular momentum gas situated in the outer regions 
of the forming galaxi es (e.g., Bek ki 1998; Rix et al. 19991 ; 
iNaab fc Burkert 200ll ; IKawata fc Gibson 2003h . On the 
other hand, the gas located at inner radii tend to sink 
toward the galactic centre where it forms new stars. T his 
outside-in formation scenario (e.g., IPipino et al. 2006bl ) is 
predicted to induce a shallow positive age gradient in the 
galactic stellar population. In Paper II we find that the 
14 low-luminosity galaxies have an average age gradient of 
0.06 ± 0.13 dex per decade. The gradients are such that 
central regions (i.e., r ^ 0.5r e ) are a few Gyr younger with 
respect to the outer parts. The stellar age of these outer 
galactic regions is on average older than 8 Gyr (see table 3 
in Paper II) In other words, the cold stellar component 
of the sample galaxies is old and it was formed at high 
redshift, i.e. z ^ 1. 

The radial variation in the ellipticity of the isophotes, 
and their prevalent discyness, observed in our data sam- 
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pie could be attributed to the cen tral concentration of stars 
produced by the sinkin g gas (e.g.. lMihos fc Hernquist 19961 ; 
iNaab fc Burkert 200 ll ). This mass congregation makes the 
gravitational potential well of the galaxy steeper, thus 
prompting a change of the stellar orbits. Specifically, it be- 
comes more axisymmetric and less triaxial. The lack of sup- 
port for box orbits makes the minor-axis tubes more popu- 
lated, thus producing discy orbits. 

The properties of the LOSVD s are also af- 
fected by the presence o f gas (e.g., INaab et al. 20061 : 
iGonzalez-Garcia et al. 20061 ). The local correlations be- 
tween /j.3,4 and v ro t/cr observed for our low- luminosity 
galaxies, is reproduced in numerical simulations by super- 
imposing a cold stellar disc compon ent in the outer region s 
(i.e., > lr e ) of the spheroidal body (JNaab fc Burkert 200lT ). 
The t ransfer of angu lar momentum between stars and gas 
(e.g., iMihos fc Hernquist 1996) and the redistribution of 
angular momentum from larger radii to inner radii produce 
a net increase of the rotational velocity of the galaxy. 

The efficiency of gas consumption is also found to 
be r esponsible in affecting the shape of the isophotes 
(e.g-. lBekki fc Shiova 19971 ). Discy isophotes are more likely 
formed in galaxies with gradual star formation, whereas 
galaxies with a rapid starburst tend to have boxy isophotes. 
The almost solar value of the [a/Fe] abundance ratios dis- 
played by our galaxies (Paper II) favours this idea. The ob- 
served abundance values are a consequence of the low star 
formation efficiency in our galaxies, thus leading to an ex- 
tended star- forming episode (^ 1 Gyr). 

In conclusion, the kinematic and photometric results of 
our analysis are in good agreement with the interpretation 
of the stellar population findings in Paper II. They suggest 
that an early star-forming collapse with a mass-dependent 
star-formation efficiency is the main mechanisms acting in 
the formation of these low-luminosity early-type galaxies. 

Galaxy mergers may also produce some of the ob- 
served feat ures. We now compare o u r results to merger sim- 
ulations bylNaab fc Burkert (20031 ), INaab et al. ( 20061) and 
ICox et al. (20061). 



Naab fc Burkert (20031 ) used a large set of collisionless 



iV-body simulations to investigate the fundamental statisti- 
cal properties of merger remnants. Simulations involve dis- 
sipationless binary mergers of disc galaxies with mass ratios 
of 1:1, 2:1, 3:1, 4:1. Projection effects are taken into account 
and the bulge components of the progenitor discs are left 
to evolve until they coalesce before the kinematic and pho- 
tometric properties of the simulated early-type galaxies are 

mea sured. 

Naab et al. (20061) re-simulated the merger remnants of 



INaab fe Burkert (20031 ) with an additional gas component 
in the progenitor disc galaxies. They considered only merg- 
ers of galaxies with mass ratios of 1:1 and 3:1. Each pro- 
genitor has 10 percent of the total mass of the stellar disc 
replaced by gas, but substantial heating processes prevent 
gas cooling and hence star formation. The gas is allowed 
to set into an equilibrium state before the disc galaxies are 

merged. 

The study of ICox et al. (20061 ) considers mergers of 
equal-mass disc galaxies. Two sets of numerical simulations 
are generated to compare the statistical properties of early- 
type galaxies formed by dissipationless and dissipational 
mergers. Progenitor galaxies can contain a purely exponen- 



tial stellar disc with a null fraction of gas, or disc where 40 
percent of the total mass is in the form of gas. Gas-rich merg- 
ers include radiative cooling of gas, star formation, feedback 
processes and induced growth of massive central black holes. 
The orbits of the merging progenito rs a re identical to thos e 
explored bv lNaab fc Burkert (200*31 ) and lNaab et al. (20061 ). 
allowing a direct comparison of results. 

In Fig. 1151 we show the distribution of the anisotropy 
parameter of our sample in comparison with the pre- 
dictions of INaab fc Burkert (20031 ). INaab et al. (20061 ) and 



ICox et al. (20061 ). The distribution of our low-luminosity 
galaxies has a median value of (i) ma i/(7o)*~ 0.86 and the val- 
ues span the range 0.4 - 2.0. The distributions of simulated 
remnants have significantly different median values for the 
degree of rotational support with respect to our galaxy sam- 
ple. However, there exists a tendency for mergers with in- 
creasingly different mass ratios to produc e more isotropically 
supp orted galaxies. The predictions of INaab et al. (20061 ) 
and ICox et al. (20061 ) show that an additional gas compo- 
nent in the progenitor discs has a significant impact on the 
galaxy dynamical support. This dissipative component in- 
crease the amount of rotation and therefore ordered mo- 
tion in the remnants. The results of a Kolmogorov— Smirnov 
(KS) test sho w that our obs erved galaxies and galaxies 
simulated by ICox et al. (20061 ) have a difference in their 
anisotropy parameter distributions significant at 4<r level. 
On the other hand, remnant s of gas-free mergers with 
mass ratios of 4:1 simulated by Naab fc Burkert (20031 ) and 
merger models of INaab et al. (2006h with mass ratios 3:1 
with progenitors having 10 percent of gas appear similarly 
distributed to our sample galaxies at almost the 60 per cent 
level. 

In Fig. 1161 the observed distribution of the isophote- 
shape parameter B4 (the discyness/boxiness parameter) 
is compared to merger predictions. The sample galax- 
ies have discy-shaped isophotes, with a distribution me- 
di an value of Ba = .30. The dissipationless models 
of INaab fc Burkert (20031 ) produce merger remnants with 
more discy isophotes when progenitors of increasingly differ- 
ent mass ratio are considere d. On the other hand, the simu- 
lations of lNaab et al. (20061 ) suggest that the presence of a 
gas component in the progenitors is responsible for changing 
the stellar orbits of the remnant. In this case, progenitors 
with a mass ratio 3:1 and 10 percent of a gas component 
produce remnants with a small amount of discy isophotes. 
A KS test marginally rejects at the 70 per cent level the 
hypothesis that dissipationless merger with mass ratios of 
2:1 and 10 per cent gas-rich mergers with mass ratios 3:1 
have the sa me distribution o f our sample galaxies. The sim- 
ulations of ICox et al. (20061 ) show that equal- mass progen- 
itors with 40 percent of gas and merger-induced star for- 
mation produce remnants with significantly discy isophotes. 
The KS test confirms that the degree of discyness of the 
simulated galaxies is inconsistent with that one of the ob- 
served galaxies. We note that the inclination angle under 
which the galaxy is observed might affect the evaluation of 
the isophotes's shape (i.e. projection effects). 

In conclusion, the kinematic and isophotal shape predic- 
tions from merger models examined are in less good agree- 
ment with our results. Although we suggest that merging is 
unlikely to be the main galaxy formation mechanism of our 
galaxies, comparison with the model predictions provides us 
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Figure 15. Frequency distributions of the anisotropy parame- 
ter for the sample of low-luminosity early-type galaxies presented 
in this work (shaded red region). The median value of the dis- 
tribution is shown as superimposed red arrow at the bottom of 
the histogram. The three series of arrows above the histograms 
indicate the peaks of the distributions for simulated galaxies . 
Bottom set of arrows: simulations from iNaab &; Burkert (20031 ) 
(labelled as NB03) for gas-free mergers with mass ratios of 1:1, 
2:1, 3:1 and 4:1. Middle set of arrows: distributions peaks from 
INaab et al. (20061 1 (labelled as NJB06) for mergers with mass ra- 
tios of 1:1 and 3:1 and progenitors with 10 percent of the total 
mass of the ste llar disc re placed by gas. Upper set of arrows: 
predictions from lCox et al. (20061) (labelled as C+06) for dissipa- 
tionlcss and dissipational equal-mass mergers; in gas-rich mergers, 
progenitors have 40 percent of the total mass of the stellar disc 
replaced by gas. 



with some constraints on the merger properties. In order to 
reproduce the observed properties the mergers are required 
to happen at high-redshift (z ^ 1), between progenitors of 
different mass ratio (at least 3:1) and with a significant frac- 
tion of the total mass in the form of gas (i.e., ^ 10 percent). 
However, this formation scenario is not able to recreate the 
the observed mass- metallicity gradient re lation of these low- 
luminosity galaxies l|Spolaor et al. 2009al ; Paper II). 

A further interpretation of our results is that low- 
luminosity galaxies were originally late-type galaxies, whose 
star formation has been truncated by removal of gas (i.e., 
strangulation) and subsequently the disc has been dynami- 
cally heated by high velocity encounters (i.e., galaxy harass- 
ment) in the cluster environment. Simulations have shown 
that late-type galaxies entering in a rich cluster can undergo 
a significant morphological transformation into spheroidals 
by en counters with brig h ter galaxies and the clust er's tidal 
field l|Moore et al. 19961 : iMastropietro et al. 20051 ). In this 
scenario, we expect the original disc to be dynamically 
heated by the interactions such that stellar orbits ac- 
quire a significant velocity component perpendicular to the 
disc. The imprint of their previous morphological nature is 
preserved in the form of an embedded stella r disc (e.g., 
|Pe Riicke et al. 20031 : IChilingarian et al. 200sl) . For exam- 



Figure 16. Frequency distributions of the coefficient B4 for the 
sample of low-luminosity early-type galaxies presented in this 
work (shaded red region). Description is as in Fig. 1151 



pie, iBeaslev et al. (20091 ) reported significant rotation at 
large radii (i.e., 4 — 7r e ) of two luminous Virgo dwarf ellipti- 
cals, using globular cluster systems as tracers of galaxy dy- 
namics. They show that the detection of such large amount 
of rotation in the outer galactic regions support the idea that 
luminous dwarf ellipt icals were originally disc galaxies. Nu- 
merical simulations IjMastropietro et al. 20051 ) predict val- 
ues of the anisotropy parameter (v/a)* similar to those of 
our galaxies, and also disky isophotes with similar B4 values. 
Moreover, the lack of counter rotation and the high incidence 
of coupled rotation between disc and bulge observed in our 
sample galaxies may favour a disc-heating scenario, whereby 
stars of the original disc contribute towards the bulge pop- 
ulation while retaining some of the angular momentum. 

However, this scenario is more difficult to reconcile with 
the stellar population properties, and in particular the mass- 
metallicity gradient relation, observed in our galaxies. We 
expect the star formation to be truncated in the disc due to 
the interactions. Thus, metallicity gradients are required to 
form in the late-type galaxies and to be somehow preserved 
during the high velocity encounters. We conclude that al- 
though the kinematic and isophotal features in our galaxy 
sample can be interpreted in the context of a morphologi- 
cal transformation from late to early types, further detailed 
numerical simulations are needed to understand if such a 
scenario can also explain the stellar population trends ob- 
served. 



8 CONCLUSIONS 

We have investigated the kinematic and photometric prop- 
erties at large galactocentric radii for a sample of 14 low- 
luminosity early-type galaxies. The radial extent considered 
in our analysis, i.e. ~ 1 — 3r e , allows us to probe more than 
50 percent of the total baryonic mass of a galaxy and more 
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than ~ 10 percent of its angular momentum. We derive ra- 
dial profiles for the kinematic parameters of rotational veloc- 
ity, velocity dispersion, and the Gauss-Hermite coefficients 
/13 and /14. We also obtain radial profiles for the photomet- 
ric parameters in both g and z-band of surface brightness, 
cllipticity, position angle and discyness/boxyness. 

We find that our low-luminosity early-type galaxies host 
an extended, cold, and old (i.e., ^ 8 Gyr) stellar component 
at radii larger than 0.5r e . This reaches out to the largest 
radii (i.e., ~ 3r e ) and the amount of rotational support 
increases as a function of the galactocentric radius. This 
implies that the low-luminosity galaxies are predominantly 
supported by rotation. The isophotes are discy-shaped and 
their degree of ellipticity is such that these galaxies can be 
well described by oblate isotropic rotator models. 

The observed properties suggest that an early star- 
forming collapse is the main mechanism acting in the for- 
mation of these galaxies. The old and extended cold stellar 
component is generated from the initial high angular mo- 
mentum gas in the outer regions of the forming galaxy. The 
low star formation efficiency of these galaxies, endorsed by 
their stellar population properties (Paper II), favours the 
observed isophotal discyness and the isophotes' ellipticity 
radial variation. Our results support and complement the 
interpretation of the stellar populations findings in Paper 
II. 

Major mergers are unlikely to be the main mechanisms 
responsible for the formation of these galaxies. In particular, 
numerical simulations of different merger categories fail to 
reproduce the amount of rotational support and isophotal 
shape observed in the low-luminosity galaxies. Although we 
can not rule out a minor merging origin for these galaxies, 
the comparison of our results with theoretical predictions 
provides us with some constraints on the possible merger 
progenitors. In order to reproduce the observed properties, 
the mergers are required to happen at high-redshift (i.e., 
z ^ 1) involving progenitors of different mass ratio (at least 
3:1) and with a significant amount of gas content (i.e., ^ 10 
percent). 

An alternative scenario is that the low-luminosity galax- 
ies were formerly late-type galaxies, whose star formation 
has been truncated by removal of gas and subsequently the 
disc has been dynamically heated by high speed encoun- 
ters in the cluster environment. However, we conclude that 
further detailed numerical simulations are needed to under- 
stand if such disc-heating scenario can also reproduce the 
stellar population properties of our sample galaxies. 
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APPENDIX A: INDIVIDUAL GALAXIES 

Al FCC 148 (NGC 1375) 

The rotation of FCC 148 increases outward to a maximum 
amplitude of ~ 89 km s _1 , where a falls to ~ 40 km s . 
The hz has the opposite sign to the rotation, often seen in 
systems with an embedded cold disk. 

The increase in rotational support outwards is also re- 
flected in the increasingly elliptical isophotes with radius, 
with the ellipticity rising from ~ 0.2 near the galaxy centre 
to > 0.6 at r > \r e (1.3 kpc). 



A2 FCC 153 (IC 1963) 

The rotation of FCC 153 increases outward and plateaus 
at ~ 133 km s , where a falls to ~ 30 km s - . The a 
profile is slightly asymmetric in the centre, being higher on 
the negative r side. The /13 has the opposite sign to the 
rotation, often seen in systems with an embedded cold disk. 
Like FCC 148 the increase in rotational support out- 
wards is also reflected in the increasingly elliptical isophotes 
with radius, with the ellipticity rising from ~ 0.2 near the 
galaxy centre to s=s 0.6 at r > 0.2r e (220 pc). The isophotes 
between 0.1 and 0.4 r e (110 pc and 440 pc) are also signifi- 
cantly disky. The surface brightness \i only slightly decreases 
with radius in the galaxy centre, and exhibits an abrupt 
break at r ~ 0.9r e , roughly at the transitional point where 
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the rotation plateaus. Thus the correlation between the sur- 
face brightness, rotation and ellipticity may be caused by a 
radially varying disk-to-bulge ratio. 



A3 FCC 170 (NGC 1381) 

The kinematics of FCC 170 exhibits several peculiarities. 
Firstly, over the monotonically outwardly increasing rota- 
tion, the rotation is enhanced in the central regions (\r\ < 
0.3 r e ), where a is also peaked, and the hz has an opposite 
sign to the rotation. This seems to suggest there is a small 
nuclear embedded disk with radius ~ 370 pc. In the outer 
parts of the galaxy the rotation is strong, at ~ 200 km s _1 , 
where the hz has an opposite sign to the rotation. This sug- 
gests substantially skewed LOSVDs at large radii, which 
could be the signature of an extended disk. This seems 
to be supported by the photometry, which shows that the 
isophotes are increasingly elliptical with increasing r. The 
ellipticity rises from ~ 0.2 at 0.1r e (120 pc) to 0.5 at 1 r e 
(1.2 kpc). 



A4 FCC 277 (NGC 1428) 

The kinematics of FCC 277 bears some resemblance to that 
of FCC 148, with an outwardly rising rotation which reaches 
~ 68 km s — 1 . There is also a weak asymmetry in the LOSVD 
indicated by anti-correlation between hz and rotation. 

The increase in rotational support outwards is also re- 
flected in the increasingly elliptical isophotes with radius, 
with the ellipticity rising from ~ 0.1 at r ~ 0.1r e (90 pc) to 
> 0.45 at r > 2r e (1.8 kpc). 



A5 FCC 301 (ESO 358-G059) 

The most notable feature of FCC 301 is an outwardly rising 
a profile and a depressed central a, which is the reverse of 
most other galaxies. The rotation seems to peak at \r\ ~ 
0.3r e with amplitude of ~ 50 km s _1 and then falls off. It is 
tempting to interpret the depressed central a as due to the 
presence of a luminous, dynamically cold sub-component of 
radius 0.3r e (289 pc) in the galaxy centre. Farther out, the 
higher a (~ 80 km s _1 ) and the lower rotation indicate an 
outward increase in dispersion support. 

The photometry supports the hypothesis of a central 
kinematic subcomponent. The surface brightness is flat or 
slightly decreasing with radius within 0.3r e (280 pc) of the 
galaxy centre, and falls off beyond that. The isophotes has a 
dramatic rise in ellipticity within the same region, reaching 
e ~ 0.6, and then falls with distance beyond that. 



A6 FCC 335 (ESO 359-G002) 

The rotation of FCC 335 is not symmetric. For positive r the 
stars' mean rotation gently rises outwards to ~ 15 km s _1 , 
while for negative r the stellar mean rotational amplitude 
reaches a maximum at ~ 30 at 0.5r e (700 pc). There is also 
asymmetry in a. In addition, the a is centrally peaked inside 
~ 0.1r e (140 kpc), where the hi is also peaked. The galaxy 
has an outwardly rising ellipticity, which reaches a maximum 
of 0.4. 



A7 VCC 575 (NGC 4318) 

The kinematics of VCC 575 is quite similar to that of FCC 
153, but with lower rotation and a. The presence of an em- 
bedded disk is also suggested 

The increase in rotational support outwards is also re- 
flected in the increasingly elliptical isophotes with radius, 
with the ellipticity rising from ~ 0.2 at r ~ 0.4r e to > 0.35 
at 0.5 > r > 2r c . The isophotes are also significantly disky 
beyond 1 r e . The surface brightness fi exhibits a break at 
r ~ 0.8r e , roughly at the transitional point where the ro- 
tation plateaus. Thus the correlation between the surface 
brightness, rotation and ellipticity may be caused by a ra- 
dially varying disk-to-bulge ratio. 



A8 VCC 828 (NGC 4387) 

The most striking feature is a distinct kinematic feature 
within 0.25 r e (200 pc) of the nucleus, indicated by a cen- 
trally lowered dispersion and an enhanced velocity gradient 
compared to the rest of the galaxy. This suggests that em- 
bedded within the slowly rotating galaxy is a bright, dynam- 
ically cold small disk component with a higher rotational 
speed. The hz does not exhibit any significant feature asso- 
ciated with the decoupled core. This may suggest that the 
central disk dominates the brightness there causing only a 
weak deviation from a gaussian LOSVD. 

The photometry also support the kinematic findings. 
The ellipticity profile also shows a break at the location 
where the dispersion changes gradient, with e roughly con- 
stant at 0.2 inside 0.2 r e , and outward rising beyond that. 
There is also a break in the B4 profile at the same lo- 
cation, indicating that the isophotes go from elliptical to 
slightly boxy from 0.2 r e outwards. Our results are in good 
agreement with the previously kinematic measurements by 
lHallidav et al. (20011 ) . 



A9 VCC 1025 (NGC 4434) 

The kinematics of VCC 1025 shows signs of an extended em- 
bedded disk, as indicated by the asymmetric LOSVD with 
opposing v and hz- At the outerparts of the galaxy v ro t and 
cr have roughly similar amplitude, indicating comparable ro- 
tational and dispersion support. Photometrically the galaxy 
is fairly round thoughout. It is possible that the disk is in- 
clined. 



A10 VCC 1146 (NGC 4458) 

The most striking feature is the presence of a compact kine- 
matically distinct core. The galaxy overall has little rotation 
and is dispersion supported. However, within ~ 0.1r e of the 
nucleus (210 pc), the rotation is significantly enhanced and 
the LOSVD is distinctly skewed, with hz having opposite 
sign to Vrot- Such feature can be attributed to the presence 
of a nuclear disk. The a is strongly peaked at the nucleus, 
and falls sharply outwards. This fall is associated with the 
sharp rise in rotation, due to the contribution of the disk 
component to the galaxy surface brightness. 

The compact kinematically distinct core is also detected 
in the photometry. There is an abrupt break in elliptic- 
ity roughly where the aforementioned kinematic transition 
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occurs. The nuclear region has a high ellipticity (e « 0.4) 
while the regions beyond 0.1r e are fairly round (e < 0.15). 
The kinematically distinct core is also slightly disky com- 
pared to the rest of the galaxy. Our results are in good 
agreement with the p reviously kinematic measurements by 
lHallidav et al. (200 ll ). 



All VCC 1178 (NGC 4464) 

VCC 1178 appears to possess a large scale disk, as indi- 
cated by a rotation curve which resembles that of a disk 
galaxy, with a sharply rising rotation which reaches a max- 
imum within 0.2r e (100 pc), and then gradually declines 
outwards. The high asymmetry in the LOSVD also supports 
this hypothesis. At r ~ 100 pc, where the rotation reaches 
its maximum, the dispersion is roughly half that of the nu- 
clear value. The kinematics is generally consistent with the 
presence of a dispersion-supported bulge plus a large scale 
disk. 

The photometry also indicates a radially changing disk- 
to-bugle ratio. The ellipticity profile has a peak at 0.2r e , 
and is lower in the galaxy center and the outer parts. In- 
terpreting the increase in ellipticity is caused an increase 
in light contribution from a disk component, both kinemat- 
ics and photometry indicate that the large scale disk has a 
maximum contribution at 0.2 r e . Our results are in good 
agreement with the p reviously kinematic measurements by 
lHallidav et al. (200 ll ). 



A12 VCC 1297 (NGC 4486B) 

The kinematics of VCC 1297 bears some resemblance to that 
of VCC 1178, with the exception that the dispersion is com- 
parably higher throughout the galaxy, the LOSVD asymme- 
try is weaker, and the decline in the dispersion profile is less 
dramatic. This may suggest a smaller contribution to the 
total surface brightness from the large-scale disk, compared 
to VCC 1178. 

The photometry shows a generally outwardly decreasing 
ellipticity with radius. The ellipticity has a "shelf" between 
roughly 0.1 and 0.4 r e , where the rotational is fastest. As 
the rotation gradually declines to an asymtoptic value, the 
galaxy becomes rounder. These seem to suggest that the 
large scale disk is more dominant at the galaxy centre than 
in the outer parts, and VCC 1297 has a higher support from 
dispersion. 



is so weak as indicated by v r ot- It is unlikely that the kine- 
matic subcomponent is a non-rotating dwarf or star cluster, 
or a disk which is orientated along the slit. It is possible 
that there exists a cold disk which is not aligned with the 
slit, such that little rotation is sampled along the slit, while 
there exists enough asymmetry in the LOSVD to give rise 
to the /13 profile. 

The photometry seems to support the above hypothesis. 
The P.A. is ~ 60 throughout the galaxy, thus indicating that 
the slit is 45 degrees offset from the major axis. The "disk" 
is associated with a region of high ellipticity (e ~ 0.5 and 
significant diskyness (B4 ~ 0.03). Farther out, the galaxy 
becomes progressively round and exhibit a high degree of 
dispersion support. 

A14 VCC 1630 (NGC 4551) 

The galaxy's rotation resembles that of a disk galaxy, with a 
maximum amplitude of ~ 70 km s . The galaxy lies close to 
the dividing line between rotational and dispersion support. 
A striking feature is the large asymmetry in the LOSVD 
within ~ 0.3r e (300 pc) of the nucleus, indicated by the hz 
which has the opposite sign to v rot . Over the same region a is 
depressed. The kinematics are consistent with the presence 
of a nuclear cold disk as well as a fainter outer disk, which 
could be the manifestation of a single structure due to a 
radially varying disk-to-bulge ratio. 

The photometry does not show any particular outstand- 
ing feature. The galaxy is moderately elliptical with elliptic- 
ity roughly fluctuating around e ~ 0.3. There is hint of a 
slightly enhanced ellipticity and diskyness within r ~ O.37-3. 
But the overall lack of a distinctive photometric feature 
seems to indicate the presence of a single embeded disk 
which give rise to the aforementioned kinematic features. 
Our results are in good agreement with the p reviously kine- 
matic measurements by lHallidav et al. (20011 ). 



A13 VCC 1475 (NGC 4515) 

The most striking kinematic feature is the lowered disper- 
sion within 0.2r e (150 pc) of the nucleus. Within this re- 
gion there is a mild counter-rotation of ~ 5 km s~ with 
respect to the rest of the galaxy, which has a maximum ro- 
tation of ~ 22 km s . There is significant asymmetry in 
the LOSVD, with hz having the opposite sign to Vrot- The 
/13 also changes sign on opposite sides of the nucleus. If the 
kinematically distinct core is caused by the light dominance 
of a kinematically cold subcomponent, then this subcompo- 
nent should exhibit quite pronounced counter-rotation. At 
first sight, it seems puzzling why the LOSVD asymmetry is 
pronounced, as indicated by the /13, while counter-rotation 



